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Abstract. Motivated by the result that an 'approximate' evaluation of the 
Jones polynomial of a braid at a 5*'* root of unity can be used to simulate 
the quantum part of any algorithm in the quantum complexity class BQP, and 
results relating BQP to the counting class GapP, we introduce a form of addi- 
tive approximation which can be used to simulate a function in BQP. We show 
that all functions in the classes #P and GapP have such an approximation 
scheme under certain natural normalisations. However we are unable to deter- 
mine whether the particular functions we are motivated by, such as the above 
evaluation of the Jones polynomial, can be approximated in this way. We close 
with some open problems motivated by this work. 



1. Introduction 

The quantum complexity class BQP consists of those decision problems that can 
be computed with bounded error, using quantum resources, in polynomial time. 
Relative to the polynomial hierarchy of classical computation, it is known that 

BPP C BQP C PP C PSPACE, 

and at the moment none of these inclusions is known to be proper [T]. Recent 
work by Freedman, Kitaev, Larson and Wang [5J has shown that the 'quantum 
part' of any quantum computation can be replaced by an approximate evaluation 
of the Jones polynomial of a related braid. A classical polynomial time algorithm 
can convert a quantum circuit for an instance of such a problem, into a braid, 
such that the probability that the output of the quantum computation is zero is 
a simple (polynomial time) function of the Jones polynomial of the braid at a S*'^ 
root of unity. For an exact statement of this see Freedman, Kitaev, Larsen and 
Wang [3], or the more detailed papers by Freedman, Kitaev and Wang [0], and 
Freedman, Larsen and Wang [71 |8J. 

It therefore follows that if we take A{L,x) to be an oracle that returns the 
evaluation of the Jones polynomial of a braid L at a point x, any BQP computation 
can be replicated by a classical polynomial time algorithm with one call to A, 
i.e. BQP C P"^. Since computing the Jones polynomial is in general a #P-hard 
problem, this does not help. However, it is not an exact evaluation of the Jones 
polynomial that is required, but an approximate evaluation at a specific point for 
braids of a specific class. Hence we may look for a weaker oracle A' such that 
BQP C P^'. 
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In a different approach Fortnow and Rogers [1] linlc quantum complexity to the 
classical complexity class GapP. In particular they show that for any quantum 
Turing machine M running in time t{n) there is a GapP function / such that for 
all inputs x 

Pr(M(a:) accepts) = 

Again evaluating a general GapP function exactly is #P-hard, however one can 
simulate M using a polynomial algorithm with access to an oracle A", where A" 
is an oracle giving an approximation to the GapP function /. 

With this motivation we examine the type of approximation needed in order to 
simulate a quantum computation, and then consider the complexity of such ap- 
proximations. It turns out that an additive approximation is sufficiently powerful. 
We should emphasize that a polynomial time additive approximation scheme is 
weaker than the familiar and much studied fully polynomial randomized approxi- 
mation scheme (FPRAS). However it is well known that any function which counts 
objects for which the corresponding decision problem is NP-complete cannot have 
an FPRAS (unless NP=RP). We show below that all functions do have poly- 
nomial time additive approximation schemes under natural normalisations. We 
also show that in two senses this is the best sort of approximation we can hope to 



achieve in polynomial time (see Theorems 4.1, 4.3 and 4.4). 



2. Quantum computing 

A link L is a smooth submanifold of S^, consisting of c(L) disjoint simple closed 
curves. A braid on m strings is constructed as follows. Take m distinct points in 
a horizontal line {pi,P2, ■ ■ ■ ,Pm) and link them to m distinct points (gi, q2, . . . , qm) 
lying on a parallel line, by m disjoint simple arcs fi in M^, so that fi starts at pi 
and ends at ^^^(j) where tt is a permutation. A braid can be closed in numerous 
ways, by identifying the points pi and qj in some way, creating a link. Similarly 
any link can be represented as a braid. In particular, the plat closure of a braid 
on 2m strings is obtained by identifying the points p2i-i and p2i, and g2i-i and q2i 
for 1 ^ i ^ m. 

2.1. Topological computing and the Jones polynomial. One of the major 
difficulties in building a quantum computer has been the sensitivity of the system 
to outside interference. Freedman, Kitaev, Larson and Wang [S] introduced the 
notion of topological quantum computing, in an attempt to make the computations 
less sensitive to small disturbances. The basic idea is as follows. One can create 
pairs of special quasi-particles, called anyons, in a 2-dimensional plane sandwiched 
between two blocks of a superconductor. The anyons have a certain probability of 
annihilating each other (leaving a vacuum) when brought together. However this 
probability changes, according to the laws of quantum mechanics when one anyon 
is moved around the other before they are brought together. Even if it is moved in 
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a complete circle around the other, on reaching its original position the probability 
of annihilation is changed. Thus a system of a large number of these particles can 
be used as a quantum computer for decision problems; pairs of anyons are created, 
moved around relative to each other, and then a predefined pair of the anyons is 
brought together. If this pair annihilate each other leaving a vacuum, this is taken 
to be an output of (or rejection), if they do not it is taken to be an output of 
1 (or acceptance). The paths in the 2-dimensional surface, combined with a time 
dimension, give rise to a 3-dimensional representation of the 'computation' as a 
braid. There remain major difficulties in constructing such a quantum computer, 
and controlling the movement of anyons. However one of the important results of 
Freedman et al. is that small changes in the paths of the anyons do not affect the 
outcome of the computation, indeed it is determined by the isotopy class of the 
braid, and therefore stable under perturbations of the paths that do not change 
the braiding itself. 

The way in which the probability changes is sufficiently subtle that such a 
quantum computer is universal in the following sense. The Kitaev-Solovay theo- 
rem [121111] together with the density theorem of Freedman, Larsen and Wang [71[H] 
yields an algorithm which given any quantum circuit on m/2 qubits and error pa- 
rameter e, outputs a braid on m strings using a polynomial number of crossings 
(polynomial in m and loge^^). The topological quantum computation using this 
braid efficiently simulates the quantum circuit, (the probability of acceptance is 
within e of the correct value). Since an algorithm for a BQP problem can be 
used to generate a quantum circuit for a given instance, the above result gives an 
explicit method for finding an equivalent topological quantum computation, and 
so the class BQP is the same under either model. 

Hence a quantum computation on m/2 qubits is approximately represented by 
a braid b. In showing that the topological quantum computation depends on the 
isotopy class of b alone [5], the following link L is considered. L is the plat closure 
of the composition of b~^,b and a small loop 7 inserted (between b~^ and b) around 
the leftmost two strings (see Fig. [ijand Fig. [2|. Both b and b~^ are needed as any 
quantum computation must be reversible; the loop 7 effects a measurement of the 
qubit represented by the leftmost pair of strings. The conclusions of [5j may then 
be summarized as the following theorem: refer to |^ for full details. 

Theorem 2.1. Let -k he a problem in BQP, with a polynomial time quantum 
algorithm A, and let J be an instance of it. For any e > 0, a link L may be 
determined in time polynomial in |J| and loge such that 

(9 9^ Pr-MH^ m 1 (, , (-l)-W+-W(-a)=^-(^)VL(e^-/^) \ 
(2.2) Pr(yi(J) = 0) - ^-p^ (^1 + ^-^^^ j <e, 

where a = e^'^l^^ and [2]^ = 2cos7r/5 and c{L),w{L) and m(L) are the number of 
components, writhe and number of minima of the link L respectively. 
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Figure 1. The braid b. 




minima 



Figure 2. The hnk L. 

The minima of [3] are the individual joins in the plat closure at the bottom of 
the braid, hence m(L) is half the number of strings in b plus one. By construction, 
the number of strings in b is twice the number of (qu)bits in the input, hence 
m{L) = |J| + 1. The writhe is defined for an oriented link, and is the number of 
'positively oriented' crossings minus the number of 'negatively oriented' crossings, 
with respect to the given orientation of L. It is easily computable. The Jones poly- 
nomial is also defined for oriented links, however the formula above is independent 
of the orientation chosen for L. Since every crossing in b appears reversed in b~^, 
these do not contribute to the writhe of L, hence the writhe is determined by the 
four crossings involving 7, and can only be —4, 0, or 4 (depending on the orienta- 
tion of the two strands passing through 7). If L* denotes L with the orientation 
of one component reversed, then VL*{t) = t~^'^^'^VL{t) , where A is the contribution 
to the writhe of L from crossings of the reversed component over (or under) the 
rest of L. Hence only reversing 7 or one of the leftmost two strings can affect the 
Jones evaluation, and it is easily checked that this is compensated by the change 
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in the terms involving the writhe. To be consistent, we will retain the notation 
[2] 5 for 2cos7r/5 from [5j throughout. 



It is Theorem 2.1 that gives rise to our interest in approximating Viit). Further 
explanation of the derivation of this equation is given in p!H] where the special but 
sufficient case w{L) = is considered (we could restrict attention to braids with 
writhe zero without affecting the main results). Although this formula involves an 
evaluation at e^'^*/^, similar results can be obtained for the n*'^ root of unity for 
any n > 5,n ^ 6 but these involve multiple L's. 

The preceding paragraphs explain how an evaluation of a Jones polynomial can 
yield the answer to a (general) quantum computation. There is a weak converse 
to this. Suppose we have a quantum computer at our disposal with which to learn 
something about a Jones evaluation of a link L. We may assume without loss 
of generality (Freedman, Larsen and Wang [7]) that our quantum computer is of 
the topological kind and thus nicely adapted to braids. We can (easily) write L 
as the plat closure of a braid b by starting with the link diagram and pulling the 
overcrossings up and the undercrossings down. Let m be the number of strands of 
this braid b. If we wish to evaluate Viia), a = e^'^*/'", we encounter an important 
constant d = 2cos7r/r. The norm |VL(a)| is bounded from above by c?™'/^ with 
|^l(«)| = d"^^'^ achieved only when L is the unlink on m/2 components, a case 
which occurs when b is the identity braid. Our quantum computer will be able to 
provide an additive approximation (see below) of |VL(a)| as a variable with range 

Given m marked points in the horizontal plane and the number a, there is a 
finite dimensional Hilbert space H on which m-strand braids act through a Jones 
representation p. The m/2 maxima (in the plat closure) determine a vector c in 
this space and the m/2 minima (in the plat closure) determine a vector in the dual 
H*, which when identified with H by the hermitian inner product, is the same c. 

We have 

^^(") /I fhW \ 
-Jr^ = (Mb)\c). 



Furthermore Prob(|0)) = |(c|p(6)|c)p, where Prob(|0)) refers to the physical prob- 
ability that below the cups, after all the 'particles' have been fused in pairs, the 
vacuum |0) is observed, that is no nontrivial particles result from these fusions. 
The last formula refiects the quantum mechanical rule that the probability of ob- 
serving an outcome, in this case |0), is proportional to the square of the component 
of the state vector in the |0)-direction. 

Because the range of |Vl(q;)| depends exponentially on the number of strings in 
the braid, m{b), our quantum computer will give much better information (sooner) 
if we succeed in displaying L with, or nearly with, the minimal m{b), called the 
braid index of L. 
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Turning to the computational question, it is a theorem of Thistlethwaite [I5] 
that when L is an ahernating hnk, with associated plane graph G, then 

VLit) = aT{G; -t, -r') 

where a is an easily computable function, and T is the Tutte polynomial of the 
planar graph. It is known p2] that even for planar graphs, computing T{G; x, y) is 
^P-hard, except when (x, y) is one of a few special points, or lies on a hyperbola 
satisfying {x — l){y — 1) = q & {I? 2}. 

Since (— e^'"*/", — e~^'^*/"), n ^ 5, is not one of these 'easy' points, exact compu- 
tation in polynomial time is not feasible (unless #P=P). It also seems unlikely 
that an FPRAS exists for these points. However the notion of an FPRAS seems 
to be much stronger than the kind of approximation that is needed in the current 
context. For any BQP language L there is a quantum Turing machine M such 
that for all x & L, M accepts with probability at least 3/4, and for all x ^ L, M 
accepts with probability at most 1/4. Therefore all that we require is to determine 
which quartile of its range ^^(e^'^*/^) lies in. We return to this topic in Section [s] 

When considering algorithms on braids or links, the size of the input is taken 
to be the number of crossings. A quantum gate on m/2 qubits is converted into 
a braid on m strings of length polylog (1/e), therefore the number of crossings in 
the braid associated with a BQP circuit is polynomially related to the number of 
gates in the BQP circuit. This in turn is bounded by a polynomial in the input 
size, hence an algorithm will either be polynomial with respect to both the number 
of crossings in the braid and the number of input qubits to the circuit, or neither. 

2.2. GapP functions. The class of counting functions which constitute #P is the 
set of functions that count certificates of membership of a language belonging to 
NP, hence #P functions are constrained to evaluate to non-negative integers. The 
class of functions GapP can be regarded as the closure of #P under subtraction, 
that is to say a function / : U h-h. Z is in GapP if and only if there exist functions 
g,h e i^P such that /(/) = g{I) - h{I) for all / G J. The class AWPP can 
be defined as follows A language L is in AWPP if and only if there exist a 
polynomial p and a GapP function g such that for all / G 3, 

4 2P(I^I) 

(2.3) I^L ^ 0«|gj5<i 

The increase in power of quantum computation over classical computation is 
that in a quantum computer there is an ability to cancel out computations paths. 
Fortnow and Rogers [1] show that this power is captured by the class GapP, 
in which a similar effect is seen. In particular they show that BQPCAWPP. It 
therefore follows that for a BQP language L, polynomial p and GapP function g 
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satisfying (2.3), determining which quartile of the range [0,2P(I^I)] contains g{I) 
would be enough to determine membership of L. 

To summarize, our foremost problems can be interpreted as finding a suitable 
approximation for the Jones polynomial of a link, Vilt), the Tutte polynomial 
of an associated planar graph, T{G]x,y), at a particular point, or for the GapP 
functions arising from BQP languages. 

3. Approximation 

Given a function : ^ for which no efficient exact evaluation algorithm 
is known, one may be interested in an 'approximate' answer instead. A stan- 
dard approach is to look for a fully polynomial randomized approximation scheme 
(FPRAS) for the problem. If ip is such a function and / G J is an input, then an 
FPRAS for ip is a randomized algorithm that given any / G J, e > will output 
e), such that 

Pr[|^(J,e)-V^(/)|>eV(/)]<l/4, 
and the running time is polynomial in |/| and e~^. 

Here one might be prompted to consider the following sort of approximation: 
suppose we know a range in which the answer lies; can we say where in that range 
the answer lies? Is it in the top or bottom half of the range, or in which quartile? 



We shall see that this approach is unlikely to be feasible, and in Section 3.1| we 
present an alternative. Clearly this type of approximation depends on the nature 
of the range. For the moment let us restrict our attention to the class of functions 
in T^P. We will make the standard assumption that for a given NDTM M there 
exists a fixed polynomial p such that for any input x, all certificates have size 
p(|a;|) (so the total number of possible certificates of M is 2^^'^'-'). We would like to 
answer the following problem, denoted by tt^: given r, for which k is the number 
of accepting certificates for x between i^^2P^I^'l) and ^2^*^l^l^? 

The problem 112 is simply to determine which inputs have more than half of all 
certificates as accepting certificates. The set of languages in this class is exactly 
the set PP of probabilistic polynomial time languages. Furthermore, 712 is clearly 
Turing reducible to Tr2s, for any positive integer s, since if n2s{x) < s then '7i2{x) = 
1, otherwise n2{x) = 2. Hence it is no surprise that this approach to approximation 
is NP-hard for #SAT, indeed the following lemma shows that any attempt to 
approximate #SAT in this way, or any problem with a parsimonious reduction to 
SAT, is unlikely to work. The proof is straightforward and we omit it, details may 
be found in [2]. 

Lemma 3.1. For A; G Z, deciding whether a CNF formula in n literals has more 
than 2""'^ solutions is NP-hard. 



When k = 1 the same decision for disjunctive normal form (DNF) formulae 
is equivalent to that for SAT, since the negation of a SAT formula is in DNF, 
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and hence for an instance F, we have #SAT(F) = 2"— #DNF(F), where n is the 
number of hterals. This observation leads to the following related lemma. Again, 
the proof is omitted and details may be found in [2]. 

Lemma 3.2. For /c G Z, deciding whether a DNF formula in n literals has at least 
2n-k solutions is NP-hard. 

This may seem more counterintuitive since not only is DNF in P, but also #DNF 
has an FPRAS [H]. On the other hand, the next lemma shows that the number 
of stable (independent) sets of vertices in a graph (#SS) can be approximated in 
this way, even though it is T^P-complete and does not admit an FPRAS unless 
NP=RP. Essentially this is because the 'natural' upper bound on the number of 
stable sets, 2", is far too big unless the graph has very few edges. For details of 
the proof see [2]. 

Lemma 3.3. Let G he a graph on n vertices. For r G Z, determining for which 
k, ^SS{G) G ^^y^2", ^2"'] is computable in time polynomial in n and r. 



Lemma |3.1 suggests that we cannot hope to fix a partition of the range and 
then determine in polynomial time in which section the answer lies; the difficulty 
associated with an NP-complete decision problem can be shifted to exactly the 
boundary between two parts of our partition. We therefore consider an alternative 
method of approximation which will meet our needs. 

3.1. Additive approximation. Our approach to approximation consists of de- 
termining a small section of the range depending on the input, and in which we 
can say the answer lies with high probability. This gives rise to an additive ap- 
proximation. 

Definition 1. Additive Approximation (AA): 

Given any function f : 3 ^ C and a normalisation u : Z+ i-^ an additive 
approximation for (/, u) is a probabilistic algorithm which given any / G J, e > 
produces an output /(/), such that 

Pr[|/(/)-/(/)l>em/|)]<l/4, 
in time polynomial in |J| and e^^. 

Note that the 1/4 in the definition could be replaced by any 6 G (0, 1/2), since 
we could reduce this error probability in polynomial time by taking several runs 
of the algorithm. Note also that most of the time we shall be considering the case 
where / is real. In contrast to the set of functions admitting an FPRAS, which is 
closed under addition but not under subtraction (e.g. #DNF(/) has an FPRAS, 
but #SAT(/) = 2" — 7^DNF(/) does not), we have the following result whose 
proof we leave to the reader. 
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Proposition 3.4. Suppose (/, u) and {g, v) admit AA algorithms, tlien tliere 
exists AA algorithms for (— /, u), {f + g,u + v) and {f — g,u + v). If, in addition, 
|/(/)| ^ and \g{I)\ ^ "^(1-^1) for all /, then there is an AA algorithm for 

The normalisation is crucial. Since we are most interested in determining where 
in the range of possible values the answer lies, we shall usually be taking u to be 
an upper bound on |/| depending only on input size. An additive approximation 
allows errors up to an absolute value of eM(|J|), whereas an FPRAS allows only 
errors up to an absolute value of e/(/). It is therefore a weaker notion of approx- 
imation, and it is easy to check that any function that admits an FPRAS also 
admits an AA algorithm under any upper bound. 

Lemma 3.5. Let f : 3 ^ be a function that admits an FPRAS, and let u : 
Z+ I— >• M satisfy |/(/)| ^ ^(l-^l) for all inputs I E 3. Then {f,u) has an AA 
algorithm. 

Note also that a given function will have an AA with respect to some normali- 
sations but not others. For example we show later that for the number of proper 
three colourings of a connected graph G on n vertices, -Pg(3) where Pq is the chro- 
matic polynomial of G, we have an AA for (Pg(3), 2*^). However for any constant 



5 > 0, (Pg(3), (2 - (5)"), does not have an AA unless NP=RP (Theorem |44|). In 
other words we can determine -Pg(3) to within an additive error e2" in polynomial 
time, but we cannot approximate to within an additive error e(2 — 5)". Note that 
if (/(/), has an AA, then for any fixed polynomial p, {f{I),u{\I\)/p{\I\)) 
also does, since we can absorb the polynomial factor in the normalisation into e 
at only a polynomial slowing of the algorithm. 

It is the determination of the 'best' normalisation for a given function that 
causes the greatest difficulties, particularly in relation to approximating ^^(t). 
Nevertheless our first positive result shows that any function belonging to ^^P 
does have an AA algorithm under very natural normalisations. 

4. Additive approximations for #P functions 

The class of functions which constitute #P can be regarded as the set of func- 
tions that count certificates of membership of a language belonging to NP. For a 
given NP-language L there will be infinitely many NDTM's which check member- 
ship of L, and the certificates for a given input I will depend on the machine used 
in verification. 

The main result of this section is that all such counting functions have addi- 
tive approximation schemes under the 'natural normalisation' associated with the 
corresponding NDTM. For example if we take f{G) to be the number of Hamil- 
tonian circuits in a graph G, then two possible NDTM's for checking membership 
of L are Mi which takes as certificates subsets of the edges, and checks that these 
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form a cycle of length \V\, and M2 which takes as certificates an ordering of the 
vertices f 1, f 2, . . . , fn and checks that the edges between any two adjacent vertices 
in the ordering, and between Vn and vi, do appear in the graph, (to avoid dou- 
ble counting we must insist that relative to some fixed ordering of the vertices 
71 : V 1, . . . ,n, we have 7r(fi) = 1 and 7i{v2) < 7r(f„)). In each case the number 
of good certificates for a given graph G is exactly the number of Hamiltonian cir- 
cuits of G, however Mi has 2^^^^^^ possible certificates, while M2 has {\V\ — l)!/2 
possible certificates. In either case the number of possible certificates is a natural 
upper bound on the number of Hamiltonian circuits. We show below that there is 
an additive approximation algorithm under the normalisation associated with any 
such bound. 

Theorem 4.1. Let f be a function in the class with an associated NDTM 
M , so that for a given instance I, M has f{I) accepting certificates, each of length 
Then there exists an additive approximation algorithm for (/, 2^*^l-^l)) that 
runs in time polynomial in \I\ and . 

Proof. Given an instance / of /, we will select t computation paths, or certificates, 
uniformly at random from the 2^^l^l) possible. We then run M using these inputs, 
and let Xi,i = 1 .. .t he indicator functions which take value 1 if and only if the 



4th 



computation path accepts I. The estimator for /(J) is then X 



Clearly E[X] = /(/). It remains to show that we can select t only polynomially 
large, such that the error bounds given in Definition [l] are satisfied. First note 
from Chebyshev's inequality that 



Pr 



|X-/(/)| ^ e2P(l^l) 



Var(X) 

e222p(|/|) 



m 

2Pi\l\) 



1- 



2P(\I\) ■ 



t 

1 

te^' 



Now if t = 4e ^ + 1, we have 

Pr 



|X-/(/)| ^ £2^(1^1)] < 1/4. 



□ 



Turning briefly to GapP functions, we get the following immediate corollary. 

Theorem 4.2. Let f be a GapP function such that f = g — h where g, h E#P. 

(1) Suppose that there are additive approximations for {g,u) and {h,v), then 
there is an additive approximation scheme for (/, maxjw, f }); 

(2) Suppose that g{I) and h{I) have certificates of length p{\I\) for all I , then 
there is an additive approximation scheme for (/, 2^). 
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Proof. From Proposition 3.4 we have that there is an additive approximation for 
(/, u+v). From Definition [l] we can halve the permitted error for only a polynomial 
increase in running time, hence there is an AA for (/, and therefore also 
(/, maxj-u, f }), which gives (i). When g and h have certificates of length p, by 



Theorem 4.1 there are AA schemes for (5^,2^) and (/i,2^), (ii) now follows from 
(i). □ 

We have seen that all functions contained in #P have an AA algorithm relative 
to normalisation by the size of the certificate space, and it is reasonable to ask if 
we could do better. However, we give two results that suggest this is already the 
best we can do in general. First we will see that sharpening our approximation to 
a logarithmic scale for the number of proper 5-colourings of a graph is NP-hard. 
Secondly, we show that the normalisation by the number of possible certificates 
cannot be improved significantly in the case of the number of fc-colourings of a 
graph. 

Theorem 4.3. Let -Pg(5) be the number of proper 5-colourings of a graph. Then 
for a general graph G on n vertices, there cannot be an additive approximation 
algorithm for (log(PG'(5) + 1), 3n) that runs in time polynomial in n and unless 
NP=RP. 

Proof. Consider a NDTM for -Pg(5) that takes as certificates any 5-co louring of the 
graph, hence the certificates are of length n [log 5] = 3n, where n is the number of 
vertices in G. We show that an AA algorithm for (log(PG'(5) + l), 3n) would be able 
to solve the NP-complete problem of determining whether a graph is 5-colourable. 

Let G be a graph on n vertices and consider the following polynomial time 
transformation. We form G~^ by adding n isolated vertices to G. If G is not 
5-colourable, then nor is G~^. However, if G is 5-colourable, each 5-colouring 
can be extended to 5" 5-colourings of G^. Therefore, if G is not 5-colourable 
log(PG+(5) + 1) = 0. Whereas, if G is 5-colourable 

log(PG+(5) + l) > log(5^5! + l) 
> 2n. 

Hence an additive approximation algorithm for (log(PG'+(5) -|- 1), 6n) could deter- 
mine whether G is 5-colourable or not in random polynomial time. □ 



Theorem 4.1 shows that for connected graphs there exists an AA algorithm 
for (Pg(/c), {k — 1)") as follows. We can take an arbitrary spanning tree on G, 
and take the set of certificates to define colourings relative to this spanning tree, 
giving k{k — 1)"~^ possible certificates. We can then adjust the normalisation by a 
constant factor (k — l)/ k. We show that this cannot be improved, in the sense that 
the normalisation (and therefore the error) cannot be reduced by any exponential 
factor. We have already noted that the normalisation can be improved by any 
fixed polynomial factor. 
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Theorem 4.4. // NP ^ RP then for any fixed /c > 3, 5 > there cannot he 
a polynomial time AA algorithm for {PG{k),(j){n)) for connected graphs G on n 
vertices, for any function (j){n) of order 0{{k — 1 — 5)"). 

Proof. Let (j){n) ^ c(/c — 1 — 5)" for sufficiently large n. Take r such that {k — 1—6) < 
{k — 1)^^^/''. Given any graph G, we form a graph H by attaching a path of length 
n(r — 1) to a vertex of G. Now 



(4.5) 
(4.6) 

(4.7) 



PG{k){k 



Pcik) 
Paik) 



PH{k) 

Pnik) 

((fc- l)l-lAOnr 

cPuik) (f){nr) 
(pinr) c{{k - l)i-iA)' 



Now suppose that there is an AA algorithm for [PH{k),(j){\H\)), then we can get 
an approximation Pnik) within an additive error of ^(j){nr). Using Pnik) and 
equation (4.6), we obtain an approximation Pcik). By equation (4.7), Pcik) is 

(j){nr) 



within an additive error of 1/2, since 

0(nr) 

C{{k - iy-i/r)nr - ^(A; - 1 - S)' 

Since Pcik) is integral, we can therefore determine it exactly. 



^ 1. 



□ 



5. Approximating VL{t) and related quantities 



We have seen in Section that our primary problem is to decide whether or 
not there exists an additive approximation scheme for (VL(e^'^*/^), [2]^^^) where L 
is the plat closure of a braid on m strings, indeed an additive approximation for 
the absolute value of the Jones polynomial suffices. We make this precise in the 
following theorem. 

Theorem 5.1. Let A he a oracle which takes as input a hraid h on m strings and 
e > 0, and returns an additive approximation for (|Vi(e^'^*/^)|, [2]^^^), where L is 
the plat closure ofh. Then BQP = P"^. 



Proof. Recall from Section [2T that given a braid on m strings, a topological quan- 



tum computer can be constructed such that the probability of output zero is 
^^^''Mm and the computer runs in time polynomial in m. Given e > 0, using 
independent runs of this computer and a standard sampling approach, the proba- 
bility of zero can be estimated to within an error of e^, where the number of runs 
is polynomial in e~^. Hence we may estimate |VL(e^'^*/^)| to within an absolute 
error of e[2]™'^^ in polynomial time. Therefore P"^ C BQP. 
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Secondly, suppose we have a BQP language and an input x. By Theorem 2.1 



we can determine a link L, of size polynomial in \x\, such that L satisfies Equa- 
tion (2.2), and the number of minima of L is |x| + 1. If x is in the language, 
Pr(0) < 1/4, hence 

(5.2) \VL{e^^'/'')\ < [2]P+'0.39. 

Whereas, if x is not in the language, Pr(0) > 3/4, hence 

1 / (_l)c(L)+M>(L)(^_^^3«;(L)y^(-g27ri/5y 



1 + [2]i V [2] 



N-i 

5 



> 3/4 



(5.3) |^L(e'"'/')| > [2]lrl+'0.65. 

Clearly use of an oracle giving an additive approximation for (|Vj;,(e^'^*''^)|, pjg^''*'^) 
will enable us to distinguish these two cases with probability at least 3/4. Hence 
BQP C P^. 

□ 

We saw in Section |2.1 the equivalence of the Jones polynomial and a specialisa- 
tion of the Tutte polynomial for alternating links, hence we would like an AA for 
a general planar graph G for 

(^(G;-e''^'/^-e-2'^'/^),M) 

where u is some reasonable upper bound. 

Hyperbolae of the form Hg := {x — l){y — 1) = q play a crucial role in the 
manipulation of the Tutte polynomial; loosely speaking the process of performing 
a tensor product on an input graph G with some other fixed graph enables us 
to 'move around' the Tutte plane. That is the new graph G ^ N satisfies: 

(5.4) T(G ® N; X, y) = f{N; x, y)T{G; X, Y) 

where / and the arguments X,Y can be computed in time polynomial in |?/| 
and |A^|. However, for any choice of A^, the new points X, Y satisfy, 

(5.5) {X-l){Y-l) = {x-l){y-l) = q. 

Thus, such transformations restrict us to remain on the initial hyperbola Hg] 
see pjOj for further details. The close relationship between points on the hyper- 
bola enables us to use an additive approximation at one point to get an additive 
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approximation at any other point {X, Y) on the same hyperbola for which there 



exists a suitable planar which transforms {x,y) to {X,Y) by (5.4). 



Proposition 5.6. Let x,y E Q and a planar graph on k vertices be fixed. 
Suppose there is an AA sceme for (T(G; x, y), u{n)) for any planar G on n vertices 
and m edges. Then there is also an AA for 



(T{G;X,Y),u{n + m{k-2)) 



where X and Y are the points determined by the transformation in (5.4) (depend- 
ing only on X, ?/ and A^). 



Proof. Let X and Y be the points satisfying (5.4). Since G ® N is planar (see the 
construction in [in|) we may use the AA scheme for 

{T{G ® N-x,y)M\V{G ® N)\)) 

to get an approximation to within an error tu{\V{G ® N)\) with probability at 
least 3/4 in polynomial time. Note that \V{G ® N)\ = n + m{k — 2). Since the 
running time of the AA scheme is polynomial in e~^, and f{N; x, y) is a constant, 
we can approximate to within an error of e\f{N; x,y)\u{n + m{k — 2)) and still 



run in polynomial time. By (5.4) we have 

T,,r, y VI T(GiaN;x,y) 
TiG.X.Y)- . 

Hence the AA for T{G ® N; x, y) yields an AA for T{G; X, Y) with error at most 
eu{n + m{k — 2)) in polynomial time. □ 

Because of the important role of these hyperbolae, it is natural to look at the 
hyperbolae containing the roots of unity (— e' "). These are Hg^, where 

g„ = 2 + 2 cos(27r/n), which cut the x-axis at 

(5.7) X = -1 - 2cos(27r/n), 

corresponding to an evaluation of the chromatic polynomial at one of the well 
known Beraha numbers Bn = 2 + 2 cos(27r/n). Since for real x and y and any 
graph N , the related points X and Y will also be real, we cannot find an such 
that we can directly relate T{G®N- 1-55,0) and T{G- -e^, -e"^). Whether 
or not we can find a point within absolute value e of (1 — -65,0) that can be 
directly related to (— e^, — e"^) is an interesting ongoing question. We present 
some positive results below, and return to these difficulties in Section [7| 



First note that by Theorems 4.1 [and 4.2 we know that T{G; x,y), a:,?/ G Z will 



have an AA scheme with respect to an appropriate normalisation, since evaluations 
at these points are GapP functions. However the drawback is that often the naive 
normalisation will be too large. This will not always be the case, for example 



APPROXIMATE COUNTING AND QUANTUM COMPUTATION 



15 



the point (1 — A,0),A G Z, gives the number of proper A colourings, and by 



Theorems |4.3| and AA here we have a best possible normalisation. 

When we consider the non-integer points, the situation is more complicated. A 
straight forward sampling approach gives the following result. 

Proposition 5.8. For rational {x, y) and a connected graph G, there exists a AA 
algorithm for the following: 

(1) {T{G;x,y),y\^\{y - l)-!^!) when {x = l,y> 1}, 

(2) {T{G;x,y),x\^\{x - when {x > l,y = 1}, 

(3) {T{G;x,y),y\^\{x-l)\^\-^) when {x>l,y> 1}. 

In other regions, in particular where there are negative terms in the expansion 
of T, cancellation between terms means that there is no longer a natural upper 
bound by which to normalise. We return to this problem in Section [7} 

6. An ALTERNATIVE APPROACH 

Returning to our original motivation, at the moment we are unable to deter- 
mine whether there is an AA algorithm for (VL(e^'^*/^), [2]^^^) where L is the plat 
closure of a braid on m strings. However we now show that in order to simulate 
a quantum computation it would be sufficient to determine the sign of the real 
part of Viie^'^^^^). Particularly in the case that the writhe of L is zero, and hence 
Viie^^^^^) is real, this appears an easier problem. 

Theorem 6.1. Let A{L) be an oracle that returns the sign of the real part of the 
Jones polynomial of the link L evaluated at e^"^^!^ . Then BQP C P"^. 



Proof. This proof follows that of Theorem 5.1 Suppose we have a BQP language 



and an input x. By Theorem 2.1 we can determine a link L, of size polynomial in 



x\, such that L satisfies Equation (2.2), and the number of minima of L is |a;| + 1. 



We now assume that w{L) = 0; the proof in the cases w{L) = 4 and w{L) = —4 



follow by a similar argument. Simplifying Equation (2.2) in the case w{L) = we 
have 

1 / (-l)^WV7(e2^*/5)\ 
(6.2) Pr(0) = -^[1+^-^ '-]. 

If X is in the language, Pr(0) < 1/4, hence 

1 ^,^(zi)!^!^^\ < 



i + [2]iV [2p-' 



l + [2] ^ 

(6.3) (-l)^(^)ri(e''^'/^) < -[2]P~'0.09 < 0. 
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Whereas, if x is not in the language, Pr(0) > 3/4, hence 



1+ ^ M > 3/4 



(_l)-(i)V^(e2-/5) > [2] 



-1 /3(1 + [2] 



2^ 



(6.4) (-l)^(^VL(e2"*/^) > [2][fl~'l.71 > 0. 

Clearly use of an oracle giving an additive approximation for the sign of Vi^e^'^^^^) 
will enable us to distinguish these two cases with probability at least 3/4. Hence 
BQP C P^. □ 

In the previous section we outlined the importance of the hyperbolae Hq : = 
{x — l){y — 1) = q to the Tutte polynomial. For x,y,X,Y and related as in 



Equation (5.4), we can determine the sign of T{G; X,Y) if we can determine the 
sign of T{G ® N; x, y). 

This gives rise to the natural question of the complexity of determining whether 
a function is greater than or less than zero, in particular the Tutte polynomial, of 
which the Jones is a specialization. It is immediate from the definitions that the 
Tutte is non-negative in the region x, ?/ ^ 0. At all other integer points on the axes 
the Tutte polynomial counts either colourings or flows, up to easy multiplicative 
factors. Since these factors may be positive or negative, we can always select one 
of either "T(G; x, y) is non-negative" or ''T{G] x, y) is non-positive" that is true, 
in polynomial time. In the above situation we are not concerned with cases in 
which the value is exactly zero, hence this would suffice. We consider the situation 
at other points in the next section. 

7. Some combinatorial and complexity questions 

We close with the following questions which have been prompted by this work. 

In Section |5] we noted that we are unable to find a suitable normalisation for 
approximating the Tutte polynomial when the expansion included negative terms. 
We return to this here and examine the chromatic polynomial to highlight the 
difficulties. We have seen that for a connected graph G, we have an additive 
approximation for (Pc'(A), (A— 1)")) for all A G However we are most interested 
in an additive approximation at the non-integral Beraha numbers. One might hope 
to achieve the above approximation for all A G M^^, however this seems unlikely 
as (A — 1)*^ is not even close to being an upper bound for Pg{^)- Indeed consider 
the complete graphs: for small 5 

(7.1) + = {l + 6){6){-l + 6)---{-n + 6 + l) 

(7.2) ^ (-l)"-2<5(^_2)!. 

This prompts the first open question. 
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Question 1. What is the best upper bound depending on A, n and m for |Pg'(A)| 
for all (planar) graphs G on n vertices and m edges? 

As far as we are aware the best upper bound known is: 

|Pg(A)| ^ |A|"-™(|A| + 1)™ AgC. 

For general graphs we can make the following small improvement. 

Proposition 7.3. Let G be a graph and let A G C, then 

|Pg(A)| ^ --1 - for-^IAI + l 

\n / \n / n 

171 

\PgW\ < (|A|)"-™(|A| + 1)" for-<|A| + l. 

n 

If G is a connected graph then 

/ m \ n-m-i / m \"^,,, , m 
\n — 1 / \n — 1/ n — 1 

TD 

|Pg(A)| ^ (|A-l|-ir— i(|A-inA| for r<|A-l|. 

n — I 

These bounds hold for all A G C, however the Beraha numbers have special 
characteristics. The evaluations of the chromatic polynomial at these points have 
some beautiful, but not totally understood, properties [16]. The values begin 
4, 0, 1, 2, 1 + r, 3, . . . and converge towards 4, where r is the golden ratio ■ The 
integers in this series are clearly central to the theory of chromatic polynomials. 
Writing P5 = 1 + r, then for any plane triangulation T on n vertices: 

(7.4) \Pt{B,)\^t'-\ 

For a connected graph G with average degree at least 3.24 (note that a planar 
triangulation has average degree 6 — 12/n), the above proposition gives. 

Hence we ask: 

Question 2. Is there a better bound for |Pg'(P„)| than there is for an evaluation 
at a general point? 

Following the results of Section |6] we are also prompted to examine the com- 
plexity of determining whether the Tutte polynomial is greater than or equal to, 
or less than zero at a given point. Recall that this decision problem is trivial for 
x,y ^ 0, and for integer points on the axes. Again considering the specialization 
to the chromatic polynomial we ask: 

Question 3. For fixed A G Q, is it NP-hard to decide whether Pg(A) is greater 
than or equal to, or less than zero? 
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Note that this is trivial for A G Z. It is also P-time decidable for A < 32/27 by 
the following theorem of Woodall [19j and Jackson [9]. 

Theorem 7.5. Let G be a graph without loops on n vertices, k components and b 
blocks. 

(1) // A < 0, then is non-zero with the sign of (—1)"; 

(2) // < A < 1, then Pg{^) is non-zero with the sign of (— l)"^*^; 

(3) // 1 < A < p, then Pci^) is non-zero with the sign o/ (— 1)"^'^^^. 

Note that PgW 7^ for A G Q\Z, since the chromatic polynomial has integer 
coefficients. It is easy to show the following. 

• Let A G Q\Z. If deciding whether PgW > is NP-hard, then it is also 
NP-hard to decide whether Pg{^ + 1) > for a general graph G. 

However since it is easy to decide for A < 32/27, the converse cannot be true for 
all A G Q\Z unless these questions are all in P. It would be interesting to know 
the answer to the following questions. 

Question 4. Does there exist a critical a > such that deciding whether PgW 
is greater than or less than zero is NP-hard for all rational A > a, A ^ Z? 

Question 5. Is this critical a equal to 32/27? 

As before we are more interested in evaluating the chromatic polynomial at the 
Beraha points than at general non-integers, and the graphs we are most interested 
in are planar. Hence we ask the specific question: 

Question 6. For planar graphs, is the problem of deciding whether Pc{Bn) is 
greater or less than zero NP-hard? 

For any graph G, not necessarily planar, it is known that Pg{B„) 7^ for 
n > 5,n ^ 6, 10, [13]. Also Tutte [16] has shown that for any planar triangulation 
the following equation holds, writing Biq = 

(7.6) Pt{B,o) = V5T'^''~'\PTmf. 

So Pt{Bio) > for all plane triangulations T, indeed a simple reverse induction 
shows that for any planar graph G, Pg(-Bio) > holds. Further, for any outerpla- 
nar graph G, we can form the planar graph by adding a new vertex adjacent 
to all original vertices. Since 

Pg+(A + 1) = (A + 1)Pg(A) 

holds for all positive integers, it holds for all A G M. Noting that Biq = P5 + 1, we 
conclude that Pg{B^) > for all outerplanar graphs G. 
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